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On the stability of hole crystals in layered cuprates. 
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Abstract. Recent STM measurements have revealed the existence of periodic charge modulations at the 
surface of certain cuprate superconductors. Here we show that the observed patterns are compatible with 
the formation of a three-dimensional crystal of doped holes, with space correlations extending between 
different Cu-O layers. This puts severe constraints on the dynamical stability of the crystallised hole 
structure, resulting in a close relationship between the periodicity of the electronic modulation and the 
interlayer distance. 

PACS. 74.72.-h Cuprate superconductors (high-Tc and insulating parent compounds) - 71.30.-|-h Metal- 
insulator transitions and other electronic transitions - 71.38.-k Polarons and electron-phonon interactions 



Since the discovery of high temperature superconduc- 
tivity, several microscopic models have been proposed in 
order to explain the complex phase diagram of the cuprates. 
In particular, the presence of an antiferromagnetic phase 
in the parent compounds has led to a huge theoretical 
effort on models with local electronic interactions, such 
as the Hubbard or t-J models. It is seldom realised that, 
when few carriers are added to an insulating system by 
doping, the long range part of the Coulomb repulsion is 
not screened, and can not be neglected a priori 

Interest in the long range interactions has been re- 
cently revived by the observation, by scanning tunneling 
microscopy, of periodic modulations of the electronic den- 
, sity of states at the surface of the cuprate compounds 



pulsion, also the interaction of the doped holes with all 
the degrees of freedom of the host material, including the 
1 — a: localised electrons in the Mott insulator, the ions 
of the host lattice, chemical impurities, etc... Instead of 
attempting this formidable task, or choosing to rely on a 
definite microscopic model, we shall approach the prob- 
lem starting from the experimental observation that the 
holes are localised in an ordered pattern, and describe the 
system in the framework of a phenomcnological Lorentz 
model. As we will show, a great insight into the prob- 
lem can be gained already at this phenomcnological level, 
leading to precise constraints on the admissible shape and 
periodicity of the charge ordering patterns. 

Let us consider the holes to be located at the sites Ri 



BiaSrsCaCuzOs-H* (Bi-2212) 



and Ca2-£!;Naa;<?frff)^l^-'^*"^ensiona/ Bravais lattice. This is a sensible 



(Na-CCOC). Although some of the features have been 
interpreted in terms of quasiparticle interference effects, 
the emergence of modulations with a dispersionless order- 
ing vector q ~ 27r/4ao directed along the Cu-O bonds, 
strongly points to the existence of an underlying charge 
order, characteristic of the less conducting (pseudogap) 
regions. Among other possibilities, the latter could be as- 
cribed to the formation of a Wigner crystal ^U] of holes 
|??1|11II12II18II14II15| — an insulating ordered state arising 
in electronic systems at low density, when the long-range 
Coulomb interactions are dominant — or to the ordering 
of hole pairs. ,16,17.iT^irn^ 

The aim of this work is to determine if the concept of 
hole crystallisation is compatible with the experimental 
observations in the cuprates. It is clear that a complete 
microscopic treatment of the problem should consider, in 
addition to the aforementioned long range Coulomb re- 



assumption because, despite the structural constraints in 
the cuprates, which confine the carrier motion into Cu-O 
layers, the long range Coulomb interactions are isotropic, 
so that the ordering within a Cu-O layer is certainly af- 
fected by the positions of the holes on different layers. 
Allowing for small (inplane) displacements Ri — > Ri -I- Uj 
around the equilibrium positions leads to the following 
quadratic Hamiltonian: 



Hr 



Y.H. + ^Y.<^ 
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El 2 2 



(1) 



The first term on the right represents all the static long- 
range Coulomb interaction effects, including the Madelung 
energy of the Bravais lattice, plus the restoring poten- 
tials acting at each Ri, due to the repulsion of the other 
holes localised at Rj ^ Hi. The second term is the usual 
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Fig. 1. The anisotropic crystal structures considered in the 
text: a) single layer compounds (d is the interlayer spacing, a 
is the hole-hole distance within a layer) ; b) bilayer compounds 
(I is the distance between neighbouring layers). 



dipole-dipolc interaction arising between particle oscilla- 



tions. The indices q;,/3 = x^y are summed, and X\ 



-a Si 






It should be noted that the terms Hi 

are screened by the static dielectric constant ts of the 
host medium, while the dipolar interactions between fast 
hole oscillations are ruled by the high frequency dielectric 
constant Coo- This discrepancy can have important conse- 
quences in the cuprates, where ej is substantially larger 
than eoo owing to the ionic polarisability. The last term 
represents the (unknown) potential which, within our phe- 
nomenological approach, accounts for the influence of the 
host material on the dynamics of the individual holes (in- 
teraction with the magnetic, ionic degrees of freedom...). 

Let us first consider a "pure" Wigner crystal, which 
is obtained by setting o^o = in the above Hamiltonian. 
PH) Among the different three-dimensional Bravais lat- 
tices compatible with the structural anisotropy imposed 
by the Cu-0 layers, we shall deliberately restrict our- 
selves to the cases which present a square charge ordering 
within the planes, as observed in experiments. Within this 
subclass, it can be demonstrated that the body centered 
tetragonal (BCT), illustrated in Fig. l.a., has the lowest 
Madelung energy. 

To determine if such anisotropic Wigner crystal is sta- 
ble, we now evaluate the frequencies WkA (-^ = branch 
index) of the corresponding collective charge density os- 
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Fig. 2. Spectrum of the collective excitations of different 
Wigner crystals, for wavevectors k within the plane. Left panel: 
the usual body centered cubic (BCC) Wigner lattice (7 = 1). 
Center panel: anisotropic BCT lattice with 7 = 0.5, which is 
unstable against shear, as signaled by the existence of purely 
imaginary solutions with u)^^ < 0. Right panel: the same 
anisotropic lattice embedded in a polarisable medium with 
rj = Eoo/es = 0.5. The eigenvalues are expressed in units of 
(jjp = A-jme^ /me^, m being the hole mass, and n the hole den- 
sity. The symmetry points indicated in the graphs are defined 
as follows, in units of n/a. Left panel: F = (0,0), H= (2,0), 



N= (1, 1). Centre and right panel: F 
X' = (1 + 72,1-72),X=(1,1). 



(0,0), N= (1 + 7^,0) 



dilations. The calculation is performed using standard 
Ewald summation techniques, at different values of the 
anisotropy ratio 7 = a/2d (a being the inplane inter- 
particle spacing, determined by the hole concentration, 
and d the interlayer distance). The excitation spectra are 
reported in figure 2. In the isotropic case (7 = 1, left 
panel), we recover the body centered cubic (BCC) struc- 
ture, which is known to be mechanically stable. Generic 
BCT crystals are also stable for weak anisotropy ratios 
within the range 0.66 < 7 < 1.07. Note that this identifies 
an interval of hole concentrations, where a square ordering 
in the planes arises naturally due to the three-dimensional 
character of the Coulomb interactions. 

For anisotropy ratios outside the range 0.66 < 7 < 
1.07, the emergence of imaginary eigenfrequencies {uj^)^ < 
0, see fig. 2, central panel) signals that the BCT struc- 
ture considered here is, in principle, mechanically unsta- 
ble, so that the preferred Wigner crystal structure has 
a different symmetry. However, energetic differences be- 
tween competing structures are expected to be very small, 
since Madelung energies are essentially determined by long 
range effects. Therefore, it is likely that the observed hole 
distributions are influenced by other microscopic mecha- 
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Fig. 3. The characteristic energy scale of the trapping po- 
tential necessary to stabilise an anisotropic hole crystal, as 
a function of the anisotropy ratio 7. Full lines corresond to 
the single layer host structure of Na-CCOC, with d — 7.75A 
and eoo = 4.5, at different values of the polarisability ratio 
(from left to right, »? = 1, 0.8, 0.5, 0.17). The dashed line is for 
the bilayer structure of Bi-2212, with Eoo = 4.9 and rj — 0.5. 
The system becomes increasingly unstable upon increasing the 
hole concentration (i.e. reducing 7) . The different periodicities 
observed in Na-CCOC and in Bi-2212 (indicated by arrows) 
correspond to a common localisation energy ~ 0.16eT/, char- 
acteristic of the host cuprate materials. 



nisms, such as the commensurability with the Cu-0 host 
lattice, which justifies our restriction to square planar or- 
derings. |23 

Let us now consider a crystal of holes embedded in a 
polarisable medium. In this case, an additional source of 
instability appears, causing an overall downward shift of 
the excitation spectrum (fig. 2, right panel). This occurs 
because, owing to the ionic polarisability (i.e. as soon as 
V = Coo/cs < 1), the static restoring potentials that lo- 
calise the holes, and their mutual dipolar interactions, re- 
sponsible for the dispersion of the collective frequencies, 
are screened by different dielectric constants. We conclude 
that some additional localising mechanism is required to 
stabilise a hole crystal in polarisable media such as the 
cuprates, and turn our attention to the full Lorentz model 
of eq. (1). The collective frequencies i7kA in this case can 
be derived straightforwardly from the solution of the pure 
Coulomb problem as: 20 



"kA 



^kA 



(2) 



Mechanical stability requires that the eigenfrequencies f2kA 
are all real. This, according to equation (2), sets a lower 
bound on the energy scale of the localising potential: luq > 
^min ; where w^j„ is defined as the modulus of the most 
negative eigenvalue uj^^ over the Brillouin zone. This quan- 
tity is reported in figure 3 as a function of 7, for different 
values of the polarisability ratio 77. We see that it is low 
and flat at large 7 (i.e. at low hole concentrations), but 
it rises sharply for 7 «< 0.5, indicating that charge pat- 
terns whose inplane periodicity is much shorter than the 
interlayer distance are very unlikely to occur, since these 
can only be stabilised by invoking an additional localisa- 



tion mechanism with a very large energy scale. This gen- 
eral observation establishes a direct link between the ad- 
missible periodicity of the charge modulations — an elec- 
tronic property — and the distance between Cu-0 layers 
— a structural constraint. To be more specific, using the 
appropriate values d — 7.75A, ao = 3.85A, Cao — 4.5 and 
77 « 0.5 [2| for Na-CCOC, it can be read directly from 
figure 2 that the observed periodicity a = Aqq (7 = 1) im- 
plies that the characteristic energy scale involved in the 
localisation of the holes is (at least) wq ~ 0.16ey. 

If the proposed picture is correct, it should apply to 
other cuprates, such as the bismuth based Bi2Sr2CaCu208+5. 
These compounds present an additional constraint related 
to the bilayer structure, for which the corresponding low- 
est energy configuration is illustrated in fig l.b. As in the 
single layer case, the calculated excitation spectrum shows 
that a pure Wigner crystal of holes is unstable due to the 
ionic polarisability, and an additional source of localisa- 
tion must be considered. Repeating the same arguments 
as in Na-CCOC, with parameters d — 15. 5A, I = SA 
and dielectric constants eoo — 4.9, 77 « 0.5,l23| we read 
from figure 3 (dashed line) that the observed periodici- 
ties a/ao = 4.3 - 4.7 pHlfelTTIIgl (7 = 0.53 - 0.58) are 
stabilised by taking cuq = 0.15 — 0.17ey, which is remark- 
ably close to the value obtained in Na-CCOC. This result 
strongly points to the existence of a common localisation 
mechanism, which could be at the origin of the different 
periodicities observed in different compounds. Following 
the same steps in the case of the single layer compound 
Bi2Sr2Cu06, which has a shorter interlayer distance than 
Na-CCOC {d ~ 12. 3A), we predict a crystal periodicity 
a/flo ~ 3.4 that would be interesting to verify experimen- 
tally 

Before enquiring about its microscopic foundations, let 
us analyse the consequences of the present scenario. First 
of all, since uq refiects the interaction of individual holes 
with the host material, in principle it should not depend 
much on the doping level. Here we take it as a constant, 
characteristic of the parent cuprate compounds. Accord- 
ing to equation (2), in the limit of vanishing hole den- 
sity, where all the lo^x —^ 0, the excitation spectrum will 
be dominated by luq. Upon increasing the hole concen- 
tration, however, the collective frequencies spread away 
from Wo due to the term w^^, which admits both positive 
and negative values. As a consequence, the longitudinal 
charge oscillations get hardened upon doping, while trans- 
verse oscillations are progressively softened, until a given 
eigenfrequency i7k=kc vanishes, leading to a polarisation 
catastrophe at a critical concentration Xc- 13,20,24,25 
The predicted softening of long wavelength transverse os- 
cillations, |13ll25j which is entirely due to the long-range 
polarisability of the medium, has been experimentally ob- 
served in systematic studies of the optical conductivity 
spectra of both electron- and hole-doped cuprates in the 
underdoped region. P^IITt] Similar signatures of the po- 
larisation catastrophe can be expected in spectroscopic 
ellipsometry and electron energy loss spectroscopy. 

In the foregoing discussion, we have assumed that all 
of the doped holes are crystallised, so that the system is 
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insulating for a; < Xc- In this case, the periodicity of the 
modulation varies continuously with x until it reaches its 
limiting value a/oo = {ni/xcY''^ at the critical concentra- 
tion {ni = 1,2 for single and bilayer host structures). [2HI 
When further holes are doped into the system beyond Xc, 
these can not be accomodated in the crystallised state, 
which is at the border of an instability. If, as indicated by 
experiments, the hole crystal survives for x > Xc, one pos- 
sibility is that the excess holes settle in the interstitials of 
the main ordering pattern, creating a superimposed mod- 
ulation, as observed in the Na-CCOC samples. |H| Another 
possibility is that the system phase separates into insulat- 
ing (crystallised) and conducting regions, as seems to be 
the case in Bi-2212.[3 The question of crystal stability in 
the presence of additional mobile charges, and the con- 
sequent screening of the long range interactions, remains 
open. 

Based on the emerging scenario, we can now specu- 
late on the microscopic mechanism underlying the locali- 
sation of the holes. First of all, the energy scale loq iden- 
tified above is definitely too large to be imputed exclu- 
sively to the binding potentials of chemical impurities, 
or to the pinning by commensurability effects. '29 On 
the other hand, the interaction with the antiferromagnetic 
background, which is ruled by typical exchange energies 
J ~ Q.leV, would constitute a viable possibility. How- 
ever, in strongly polarisable materials, the excess charges 
added by doping are expected to form dielectric polarons. 
Indeed, it can be recognised that the value of wq derived 
in this work coincides with the locus of a well defined 
absorption band, ubiquitous in the infrared optical spec- 
tra cr(w) of strongly underdoped cuprates, which is gener- 
ally ascribed to the formation of polarons, [3D] suggesting 
that the additional mechanism required to stabilise the 
hole crystal is the polaron self-trapping votentia,l.\^M\ The 
Wigner crystallization of polarons in the insulating phase 
of the cuprates has been proposed independently by Re- 
mova et al.,PJ| and Quemerais. ^2] The quantum melting 
of such polaron crystal has been studied in refs. |13II25II32[ 
I33j yielding similar conclusions as in the present paper. 

Our analysis demonstrates that the periodic modula- 
tions observed at the surface of the cuprates are compati- 
ble with the crystallisation of holes, arising due to the com- 
bined effects of the long range Coulomb interactions, and 
of some additional localising phenomenon, whose charac- 
teristic energy scale is ujq ~ O.WeV. Although several mi- 
croscopic mechanisms can be involved, polaron formation 
appears as a good candidate to stabilise a hole crystal in 
the cuprates. 

The authors thank D. Mayou and B.K. Chakraverty 
for useful discussions. G.R. acknowledges kind hospitality 
at LEPES, Grenoble. 
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